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AN INTERESTING LISP FUNCTION 


Ikuo Takeuchi (1978) of the Electrical Communication Laboratory of Nippon 
Telephone and Telegraph Co. (Japan's Bell Labs) devised a recursive function 
program for comparing the speeds of LISP systems. It can be made to run a long 
time without generating large niimabere or using much stack. The program is 


_tak(z, y, z)< «— if z <y then y 
else taktak — 1, y, z), tak(y — 1,z, z), tak(z — 1, x, y)), 


_ where the variables may range over the integers (including negative) or else over 
real numbers. The program has similar properties in the two cases, but proving 
termination seems more tedious if arbitrary real numbers are allowed. The program 
has several interesting features. . 

1. Inspection suggests that tak satisfies the equation 


tak(z +a, y + a,z-+ a) = tak(z, y,z) + a, 


whenever the computation terminates, and this can be shown by subgoal induc- 
tion. Namely, it is true for the non-recursive case, and assuming it for the referred 
sets of arguments yields it for the main set. A formal al can proceed along the 
lines suggested in (McCarthy 1978). 

2. Experiment using LISP indicates first of all that the value of tak(z, y, z) is 
always one of z, y or z. I don’t see a proof of this fact in isolation. 

3. Like the 91-function, the program computes a simple non-recursive func- 

_ tion. Using LISP to compute some values of tak leads to the guess that it is the 

same as l 


tak0(z, y, z) = ifz < y then y else if y < z then z else z. 


Substitution shows that tak0 satisfies the functional equation for tak, and therefore 
by the minimization schema of (McCarthy 1978), 


tak(x, y, z) = tak0(z, y, z) 
whenever the former terminates. A fortiori, this establishes that tak(z, y, z) takes 
one of the variables as value, but maybe that fact could be proved separately. 


4. In order to prove that tak is total, we write a “derived program” dtak(z, y, z) 
that computes the depth of recursion involved in computing tak(z, y, z) using call- 


http:/Awww.artinfo-musinfo.org Lisp Bulletin #3, December 1979, page 8 / 55 


7- 


by-value. We have 


=“ yz) + if z < y then 0 
else 1 -+ max(dtak(z — 1, y,2), 
dtak(y — 1,2, 2), 
dtak(z — 1,2, y), 
dtak(tak(z — 1, y, z), tak(y — 1,2, x), tak(z — 1,z,y))). 


Experiment with LISP leads to the conjecture that for integer values of the vari- 
ables, dtak is extensionally equivalent to 


dtak0(z, y,z) = dtak00(c — y,z—y), 


where 


dtak00(m,n) = if m <0 then 0 
else if n > 2 then m + ne 1)/2—1 
else if n > 0 then m 
else if n = —1 then (m + 1)(m + 2)/2—1 
else (m — n)(m— n+ 1)/2— m — 1. 


We don’t bother to verify the conjecture. Instead we use dtak0 as a rank function 
to prove Vi.®(i) by course-of-values induction where 


(i) = Vryz.(dtak0(a, y,z) = i D tak(z, y, z) = tak0(z, y, z)). 


Since we already know that tak0 satisfies the functional equation of tak, we need 
only show that in the recursive case of tak, i.e. when z > y, the referred arguments 
pi of lower rank than the main ones. Thus we must show that each of dtak0(z — 

1, y,2), dtakO(y — 1,2, 2), dtak0(z — 1,2, y), and déak0(tak0(¢ — 1, y, 2), takO0(y — 
1,z, x), takO(z — 1,2, y) is strictly less than dtak0(z, y,z). Each inequality follows 
from a separate easy case analysis. Presumably termination could be proved for 
the non-integer case by a similar argument with a more complicated formula for 
diak00. We leave this as an exercise for the reader. 

5. Like Morris’s program 


morris(z, y) — if z= 0 then 1 else morris(x — 1, morris(z, y)), 
tak is more quickly computed by call-by-need (Vuillemin 1973). In fact the num- 
ber of function calls appears to be exponential with call-by-value and quadratic 


with call-by-need. The culprit is the third argument of the outer call, namely 
tak(z — 1,«, y) which is often unneeded. Unlike morris, however, tak is total. 
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A call-by-need version of tak is given by 
ntak(z, y,2) — vtak(z, y,2), 


where 


viaku+— if numberpu then u 
else if n d u then viak( a u)— 1 
else {vtak a u, vtak ad u}[Azy. 
if z < y then y 
else vtak((z — 1, y, add u}, (y — 1, add u, 2), ( add u— 1, x, y))). 


vtak is obtained from tak in a somewhat ad hoc way. Since we don't always 
compute all arguments of a function we must work with triples whose elements are 
either numbers or applications of functions to arguments. The only functions that 
occur are vtak itself and subtracting one. Therefore, a number stands for itself, 
an application of vtak is represented by a triple, and an application of subtracting 
one is represented by a list of one element—the inner expression from which one 
is to be subtracted. Perhaps I’m being thick, but it seems to me that call-by-need 
requires lists or at least putting symbols on the stack. 

MACLISP versions of tak and friends are in the file TAK.LSP[F78,JMC] at 
SU-AI. I thank Don Knuth for suggesting call-by-need. The external or publica- 
tion LISP notation is as in (McCarthy and Talcott 1978). It has the following 
features: (1) a and d are used for car and cdr, and their compounds are formed 
similarly. (2) The infix . is used for ‘cons, and (T, y,...,2) is used for list[z, y,..., 2]. 
(3) {z, y,...,2z}f is used for f[x, y,...,2] whenever we prefer to write the arguments 
first and the function name later. 
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